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Abstract: We provide an exact analytical calculation at the two-loop level in the 
abelian limit of the leading power correction for the C parameter in e+e~ annihi- 
lation. We compare our results to the numerical value obtained employing the soft 
approximation, the abelian part of the Milan factor. We demonstrate that a simple 
proportionality holds between the leading power corrections to the C parameter and 
to the longitudinal cross section in the soft region, and we verify that this proportion- 
ality holds for the full two-loop abelian contribution computed here. We comment 
on the possibility of extending this technique to other event shape variables and 
distributions, as well as to the non-abelian contributions. 
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1. Introduction 

The study of event shape variables in both e^e~ annihilation and DIS is at present one 
of the most active areas of phenomenological investigation in QCD. From a theoretical 
viewpoint existing predictions for the distributions of these variables emerge from the 
application of several of the most advanced techniques of perturbative QCD, hence it 
is particularly satisfying to be able to confront theory successfully with experiment. 
Theoretical predictions for power corrections to shape variable distributions in fact 
require all order resummations of perturbation theory, either within the framework 
of soft gluon resummation [1-5], or in the context of renormalon techniques [6- 
10], extended to account for the necessary two-loop effects [11-14]. In either case, 
the theoretical framework is adjusted for the inclusion of non-perturbative power 
corrections, using methods and models that may ultimately pave the way for a better 
understanding of low energy behaviour of QCD [5,10]. On the experimental side, 
these variables have proved relatively simple to measure and have been one of the 
most popular sources for the extraction of the strong coupling from fits to the data 



In the present paper we will make use of the dispersive approach developed 
in [6-10]. In its original form, the dispersive approach was based on the resumma- 
tion of fermionic bubble chains, a procedure which is strictly consistent only in the 



[15-18]. 
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abelian limit, but can be reformulated under rather general assumptions in terms 
of an effective low-energy QCD coupling, for which a dispersive representation is 
adopted. Since the dispersive variable plays a role which is formally similar to a 
gluon mass, this approach is sometimes referred to as the "massive gluon scheme" . 
It was soon realized (first in a study by Nason and Seymour [19], who considered 
the thrust variable) that this method could not directly be applied to general event 
shape variables, which arc sensitive to the details of the decay process of the par- 
ticular gluon being treated as "massive". A method to deal with the problem was 
proposed in Refs. [11,20], where it was shown that the bubble resummation could 
be performed without integrating inclusively over the gluon decay products, and a 
calculation of the leading power corrections to the longitudinal cross section was 
performed in the abelian limit. A much more general analysis of the effects of non- 
inclusiveness on the most commonly used event shape variables was performed in 
Refs. [12, 13], where the two-loop correction due to gluon splitting was computed 
including non-abelian terms. The main result of this detailed analysis is that the 
whole effect of non-inclusiveness, within the framework of the dispersive approach, is 
largely universal for the shape variables studied. Each observable acquires a two-loop 
enhancement factor to the naive "massive gluon" calculation, and this enhancement 
factor, now called the Milan factor, is observable independent. All dependence on the 
chosen observable is encoded in the massive gluon result, so that the phenomenol- 
ogy is basically unaffected. Moreover, shape variables in DIS receive an identical 
enhancement [14]. 

The calculations leading to the Milan factor are done using soft two-loop matrix 
elements. This is justified because shape variables, being by construction essentially 
linear in particle momenta, are expected to receive leading (1/Q) power corrections 
only from soft gluon emission. In addition, soft kinematics is employed, where ef- 
fects such as terms involving the square of the small transverse momentum of the 
gluon decay products are discarded, since they can only contribute at the level of 
l/Q^ corrections. The universality of soft gluon radiation coupled with an underly- 
ing geometrical universality (linearity in emitted transverse momenta) in the shape 
variables themselves thus leads to the universal Milan factor appearing in every case. 

In the present paper we test some of these assumptions by considering the 1/Q 
correction to two event shape variables, the C parameter and the longitudinal cross 
section. In the case of the C parameter we perform the full calculation in the abelian 
(large n/) limit, in a manner analogous to the longitudinal cross section calculation 
of [11], and we arrive at an analytical result for the 1/Q behaviour. The correction 
to the C parameter was calculated in the soft approximation in Ref. [13], and the 
Milan factor enhancement obtained for it. The results we get here should be directly 
comparable to the abelian part of the Milan factor, since we allow only for gluon 
splitting into quarks. 

We also show that a simple relationship holds between the C parameter and the 
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longitudinal cross section in the soft approximation. This relation is respected by 
the 1/Q corrections computed here for C and in reference [11] for the longitudinal 
cross section, suggesting that the 1/Q corrections should be correctly obtained from 
the simplified soft approximation advocated in the computation of the Milan factor. 

We find however that the numerical result for the two-loop enhancement factor 
for the C-parameter in the abelian limit disagrees with the abelian limit of the 
Milan factor: if we take the ratio of our result to the massive gluon one, we find an 
enhancement factor of 157r^/128 = 1.157, whereas the corresponding enhancement 
factor computed numerically in Ref. [13] is of the form 1 + r'^ni — 1-078, where rfj""* is 
the "non-inclusive" piece of the Milan factor in the abelian limit. ^ In essence, the 
two results would agree if the correction factor r[^^ were doubled. 

The layout of the present paper is as follows. In Section 2 we introduce the 
observables C and (7^, and we demonstrate their equivalence in the soft limit. In 
Section 3 we show our calculation of the C parameter in some detail. In the final 
section we compare our results to those obtained by other authors and briefiy discuss 
the different methods, as well as possible developments and applications. 

2. Definitions and the soft approximation 

2.1 C parameter 

The C parameter [21] can be defined as 

where is the photon four-momentum, = Q"^, and the sum is over all outgoing 
particles, so that each pair of particles is counted twice. An equivalent definition, 
valid in the e+e~ centre of mass frame, is 

C = l±\PML^n-9.^ , (2.2) 
i,j=i ^ 

where the sum once again runs over all outgoing particles. From these definitions, 
it is clear that the C parameter vanishes at the Born level. The first contribution 
within perturbation theory appears at 0{as) when one has a three body final state. 

2.2 Longitudinal cross section 

The longitudinal cross section can be defined, within the context of single particle 
inclusive annihilation, by considering the differential cross section for the production 



^To arrive at this number for r^"^ we take Ca = in Eq. (4.16) of Ref. [12], in the numerator as 
well as in (3o. 
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of a hadron with given energy and angle with respect to the beam axis. The angular 
dependence of this cross section can be organized in the form [22] 



dxd cos 6 



- ( 1 + cos^ e ] + - sm^ ^— ^ + - cos ^ 
8 ^ dx 4 dx 4 



da\ 



dx 



(2.3) 



with X — 2ph ■ q/q^ being the hadron energy fraction and 9 its angle with respect to 
the beam direction. The three terms on the right-hand side are referred to as the 
transverse, longitudinal and asymmetric contributions. The first two contributions 
respectively arise from gauge boson polarisation states transverse and longitudinal to 
the direction of the outgoing hadron h; the last term is a parity violating contribution 
that is absent in purely electromagnetic annihilation. 

The longitudinal cross section itself is usually defined as the first moment of 
da^/dx, and it can be projected out from the single particle inclusive differential cross 
section in Eq. (2.3) by multiplying with a suitable weight function and integrating 
over cos 9 and x, according to 



Note that since we are measuring an inclusive cross-section we can replace the 
hadronic sum above by the corresponding partonic one. 

In the next subsection we shall make use of the above definitions to write down 
expressions for the longitudinal cross section and the C parameter in the soft ap- 
proximation, relevant to the computation of l/Q corrections. 

2.3 Soft approximation 

At the partonic level, in the soft approximation, the annihilation process produces 
a pair of hard fermions ("primary" quark and antiquark), dressed by soft gluon 
radiation. The primary quark and antiquark arc essentially back to back and have 
an angular distribution about the beam axis that may be approximated by the purely 
transverse 1 + cos^^^o pattern (the first term in Eq. (2.3)), where 6q is the quark angle 
with respect to the beam. In this approximation one is neglecting the recoil of the 
hard fermions due to soft gluon emission. 

Let us consider the contribution to ai of a soft gluon, which is known to be the 
term responsible for the appearance of a l/Q correction [11,23]. In order to do this, 
we need to introduce the gluon angle with respect to the beam 9g, while the gluon 
direction with respect to the parent quark may be specified by a polar angle 9 and 
an azimuthal angle (p. ^ In terms of these variables, the soft gluon contribution to 

^An identical calculation relevant to the flux tube model of hadronisation was described by 
Nason and Webber [24] . 




dxd cos 9 



(2.4) 
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the longitudinal cross section can be written as 



1 /"i /"i /■^'^ d(f) / \ 
= - dxx / d cos 0„ d cos On d cos Of — (2 — 5 cos^ Oq ) 

2 Jo i-i Jo 27r ^ ^ 

3 d^a^ 

X -(1 + cos^^o)^ — ; -5 (cos^p — cos^o cos^ + sin^osin^cos^) (2.5) 

8 dxd cos 6 

— - [ dxx [ 

4 Jo J-i 



acos^sm B- 



dxd cos 9 



where x is the gluon energy fraction, x = 2Eg/Q, while here d^cr^ /dxd cos is the 
differential cross section for soft gluon emission, given in the soft approximation by 
the standard qq antenna pattern. 

The relationship with the C parameter is easily established by noting that the 
only terms in the sum over outgoing particles that will contribute to the 1/Q cor- 
rection for the mean value of C are the terms pairing one of the primary fermions 
with the soft gluon. The contribution from the pairing of the primary quark to the 
antiquark is so that it can be discarded in the back-to-back (soft) limit. One 

finds then that in the soft approximation, concentrating on the terms responsible 
for 1/Q corrections, the mean value of C can be calculated by weighting the matrix 
element d^a^ /dxd cos 9 with the function 

C('/^^^^^^E,sin'9^^xsm'9 , (2.6) 

which is precisely six times the weighting function for ctl as given in Eq. (2.5). One 
concludes that the power corrections to the expectation value of C and to 0"^ are 
proportional, according to 

^C)WQ) ^ Q a^/"^) . (2.7) 

In the next section we will calculate the exact value of the power correction for 
the C parameter, taking proper account of gluon decay but in a large Uf limit. We 
will then be able to compare with the calculation of [11], and verify that the 1/Q 
corrections to (C) and ai indeed obey Eq. (2.7). 



3. Calculations 

In this section we provide some details of our calculation of the average (C) due to 
gluon splitting into a ( "secondary" ) quark-antiquark pair. We will begin by perform- 
ing the calculation in the naive massive gluon scheme, although such a calculation 
is manifestly incomplete, since we shall need it in order to compare our results with 
other available calculations. 
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3.1 Massive gluon scheme 

In the dispersive approach (for a review, see [25]), for a sufficiently inclusive ob- 
servable, defined as an observable insensitive to the momentum distribution of the 
gluon decay products, it is possible to classify power corrections by computing the 
observable at one loop with Feynman rules appropriate to a massive gluon. One 
can then show [6-10] that power corrections are in one-to-one correspondence with 
nonanalytic terms in the expansion of the result in powers of the squared gluon mass. 
It is easy to see that the C parameter, like most event shape variables, is sensitive to 
the details of gluon decay, and thus does not belong to the class of observables that 
can be dealt with using this method. One can however use the massive gluon result 
as a first estimate of the size of the correction, and parametrize the full result as 
a (perturbativc) enhancement of the massive gluon calculation. This is particularly 
useful since one can argue [12, 13] that all dependence on the particular observable 
at hand is contained in the massive gluon result. 

Let us then consider a. q,q,g final state and let us label the four-momenta of the 
quark, antiquark and gluon by pi, p2 and k respectively. The C parameter assumes 
the simple form 

^^^(1-X0(1-X.)(1-X3) 
X1X2X3 

where xi^2 — 2pi,2 ■ q/q^ are the energy fractions carried by the quark and the anti- 
quark respectively, while x^ — 2k- q/q^ is the gluon energy fraction. Energy conserva- 
tion implies J^i^i = 2. Notice that although we are doing a massive gluon calculation 
we have discarded the gluon mass in the definition of the event shape. This approx- 
imation has no effect on the coefficient of the 1/Q correction, in the present case; 
moreover, in order to compare our results with the Milan factor, we need to adopt 
the massless definition of C in accordance with the procedure of Refs. [12, 13]. 
The expectation value of C in the present approximation is given by 



{C)(mg) = — [ da C 

(To J 



where ctq is the Born cross section, e is the square of the gluon mass divided by g^, 
and {—ie){—ig)W^" is the matrix element for the decay of the virtual photon (with 
polarisation index ^) into a quark, an antiquark, and a gluon with polarisation index 
a. Current conservation then implies that g^VT^" = and kaW^"" = 0. In terms of 
the energy fractions Xi, 



xl + xl + 2e(xi +X2 + e) 



{i-x,){i-x2) {i-x,y {i-x2y 



. (3.3) 
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The integrals in Eq. (3.2) are easily computed. Expanding the result in powers 
of e we get 

(C)(mg) = ^«s[47r2- 33 - 127rV^+0(e)] . (3.4) 

In the language of the dispersive approach, this translates into a power correction, 
generated by the term, given by (see [10]) 



{C)IZ = ^^^ , (3-5) 



where the dimensional parameter Ai is interpreted as the first (half integer) moment 
of the non-perturbative component of the effective coupling, as defined in [10], 

Ai-^r^P^ SaAl^') ■ (3.6) 

A completely equivalent expression for the power correction follows from a purely 
perturbativc analysis, as a consequence of the sum over "bubble" graphs (sec [6,7]). 
In that case the coefficient Ai is related to a slightly different effective coupling 
function, related to the present one essentially by a change of renormalization scheme 
(for a comparison of the two points of view, see [25]). 

3.2 C peirameter with gluon splitting 

Let us now consider the splitting of the gluon into a quark-antiquark pair in more 
detail. We focus on a four-particle final state with a "primary" qq pair, carrying 
momenta pi and p2, while the secondary quark and antiquark have momenta ki and 
k2. Summing the relevant graphs, and summing in each graph over all insertions of 
fermion bubbles in the gluon propagator, we are led to replace Eq. (3.2) by 

X j rfLips[A; ^ fci, fca] Tr[7" /(fi7^ ^2] C . (3.7) 

Here (iLips[q' — > {pi}] is the appropriate Lorentz-invariant phase space measure, 
n(/c^) is the renormalized one-loop gluon vacuum polarization induced by quarks, 
and L^'^ is the leptonic tensor, which will be substituted by the corresponding average 
over the beam orientation, (L^*^) = —4Q'^g^^'^/3, as customary when working with 
event shape variables. In Eq. (3.7) we have factorized the four-particle phase space 
to introduce an explicit integral over the square of the gluon four-momentum, k"^, 
which of course plays the role of the dispersive variable ( "gluon mass" ) in the present 
calculation. Thus we will be interested in the /c^-dependent integrand of Eq. (3.7), 
and in particular in its expansion in powers of V^- Notice that a factor of 1/(1 + 
Il{k'^)) is just what is needed to turn the (fixed) coupling ag into the running coupling 
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evaluated at scale k^, in the abelian limit. For C we take the full expression for 4 
outgoing particles, which we write as 

C = C(p) + C(^rn) + C{s) : (3.8) 

where C(p) is the "inclusive" term, involving only the momenta of the primary 
fermions, 

C(s) is the term involving only the momenta of the secondary fermions, 

and finally C(m) is the sum of all mixed terms. Using the symmetries of the integral 
in Eq. (3.7) we can freely replace k2 with ki and p2 with pi in C(^rn), thus we can use 

Since we are interested in the distribution of C as a function of the gluon "mass" k"^, 
for each of the three contributions to (C) we write 

/C)u)^ f— ^ , (3.12) 

^ J e |l + n(g2e)|2 WV J , y J 

where i — p,m,s and e = kP'jCf'. Our task is to compute the distributions C(i)(e), 
and extract the non-analytic behaviour for small values of e. 

The primary contribution to C is clearly the easiest to evaluate, since the in- 
tegration over kx and /c2 can be done inclusively. Using the transvcrsality of the 
hadronic tensor, kaW^"" = 0, one can simply substitute for the bubble integral 

J dLips[A; ki, /cg] Tr^ja hip ¥2] Q^(^9ap ■ (3.13) 

One finds then 

The remaining integrals are easily performed, and expanding around e = yields 



C(p)(e) 



CpaiTRnf 
27r 3n 



1 "^"^ 

-Sine- — + e»(e) 



(3.15) 



Notice that C(p)(e) by itself is infrared divergent. The divergence will be cancelled 
by the mixed contribution, to which we now turn. 
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The e-distribution arising from the mixed terms is by far the most difficult to 
evaluate analytically. It is given by 

3 alTnUf 1 n-e ri- 



Nc 2% Q^e Jo 



dxi 



l—xi 

l—x\—e 



1 



dx2 -^W^^WfT^f, 



Pi ■ Q 



where 



dUps[k ki, k2] Tr[7„ /^i7/3 



(3.16) 



(3.17) 



There are several ways to evaluate the bubble integral T^/s. Perhaps the most 
straightforward is to note that Tap is symmetric, and obeys k'^T^p — 0. Then we can 
write the decomposition 



a/3 



^1 Q \9al3 TT- + ^2 Pla T^^^a [ Pip —kf. 



+ A3 



Pla 



k^ 

Pi ■ k 
k 



2 



A;2 

Qp - ^-j^^p] + (a ^ /5) 



A;2 



A ( Q- k, 

+ ^4 Ua - -p-^c 



IP 



q ■ k 



kf: 



(3.18) 



where the Ai are scalar functions of pi, A; and q, which can be evaluated, for example, 
by integrating Eq. (3.17) component by component in the rest frame of k. 

It follows from Eqs. (3.16) and (3.18), and from current conservation, that only 
four independent projections of the tensor W^°'W*^ are needed. The projection with 
the metric tensor g^p gives Eq. (3.3); the other projections are given by 



PiaPipW^^'^W*/ = AN.CfQ 
{piaQp + qaPip)W'^''W*/ = 8N,CfQ^ 

+ 



2x j{x3 - 1 - e) 

xjjxs - 1 - e) 
(l-xi)2 

xKxo, — 1 — e) 
xlix3 - 1 - e) 



+ 



(3:3-1-6)^ " 
(1 - a;i)(l - X2)_ 

(l-a:i)(l-a:2) 



(1-X2)2 

To simplify the computation, a convenient change of variables is 



u — Xi — X2 

e ^ 
V -\ — — 2 — xi — X2 

V 



X3 



(3.19) 



(3.20) 
(3.21) 



In terms of these variables, the scalars entering the decomposition of the tensor Tap, 
Eq. (3.18), have the general form 



Q"^ . _ Pi{u,v,€) + Qi{u,v,e)lav + Ri{u,v,e)\ne 



Pi ■ Q 



(v^ — e)^(v^ — 2v — uv + e) 



(3.22) 
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where Pi, Qi and Ri are polynomial functions of u, v and e, whereas the projections 
of the hadronic tensor given in Eq. (3.19) are rational functions of the same vari- 
ables. All the integrals involved in the computation of Eq. (3.16) can be performed 
analytically, through a rather lengthy process of partial fractioning. The result is 
expressible in terms of di- and trilogarithms of functions of e, and can be expanded 
around e = yielding 



Sine + 47r2 - - - + 0(e) 



(3.23) 



As promised, the infrared divergence (as e — > 0) of Eq. (3.23) cancels the one of 
Eq. (3.15). 

The final contribution to (C) is that from the secondary terms, arising from the 
distribution C(s)(e). By inspection, this contribution must be at least 0{e), since 
the corresponding event shape, Eq. (3.10) is quadratic in the gluon energy in the 
soft region, and thus cannot contribute to the leading power correction. This can be 
confirmed by an explicit calculation along the lines of the ones leading to Eqs. (3.15) 
and (3.23). 

Our final result for the C parameter e distribution is 



_ Cf alTRUf 
~ 27r 37r 



47r^-33-^v^+0(e 



(3.24) 



Following [6,7,11], one can directly compare Eq. (3.24) with Eq. (3.4), upon noticing 
that the coefficient a'^Tjin f / (Sn) should be read as a^/Sg, where Pq is the fermion 
contribution to the one-loop (3 function. From here the power correction can be 
obtained by using the renormalisation group equation for the running coupling to 
replace the "spectral function", PQa^^k"^), by the logarithmic derivative of the cou- 
pling with respect to its argument and integrating by parts. The coefficient of ^/e is 
mapped to the coefficient of the 1/Q correction as in Eqs. (3.4) and (3.5). 

Using the language of Refs. [12, 13], one would interpret Eq. (3.24) as the sum 
of the "naive" contribution to the series of power corrections plus the fermionic part 
of the "non inclusive" correction to it. We find that the exact calculation leading to 
Eq. (3.24) gives a 1/Q correction larger than the one computed in the massive gluon 
scheme by a factor 157r^/128 = 1.157, as announced in the introduction. 

3.3 Comparison with 

A calculation analogous to the one just outlined was performed for the longitudinal 
cross section in Ref. [11], where a model for the distribution in energy fraction of 
the power correction was also constructed. It was also noted there that the 1/Q 
correction to arises from a resummation of power corrections of the form l/{xQy^ 
in the distributions, which have no correction with odd powers of 1/Q; this was also 
the conclusion reached in [23], within the massive gluon model. This shows that 
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the leading power correction to the total ai is entirely due to the soft gluon region, 
a; — > 0. The result given in Ref. [11] for ai can be written in the form 

IStt^ 



(7L(e) 



CpaiTjinf 
2tv 3tt 



1 - 



64 



(3.25) 



One sees that the two results indeed satisfy Eq. (2.7), which was obtained in the soft 
approximation. This confirms that the leading power correction arises entirely from 
the emission of soft gluons, and from their subsequent splitting. 



4. Discussion 

We performed an analytical calculation, with two-loop accuracy and in the abelian 
limit, of the leading power correction to the C parameter measured in e'^e" annihila- 
tion; we showed that our result is simply related to the corresponding result for the 
longitudinal cross section, and we compared it to existing calcTilations. We explicitly 
checked that the entire power correction comes from the region in which the splitting 
gluon is soft: this can be shown by tracking the contributions to the final answer 
through the calculation, and is confirmed by the fact that the simple relationship we 
find between (C) and is a property of the soft approximation. Our calculation 
is thus a positive test of the applicability of the soft approximation to the compu- 
tation of 1/(5 corrections to event shape observables. We find however a numerical 
discrepancy with the results of [12,13], which we should reproduce in the abelian 
limit. 

Our technique, in essence a straightforward if lengthy two-loop calculation, is 
applicable in principle to all event shape variables. Although in some cases it may 
turn out to be too cumbersome to generate the full analytical result, it is always 
possible, and in fact simple, to produce a one-dimensional integral representation 
for the answer, from which the coefficient of the desired power correction can always 
be derived numerically. The fact that the calculation is performed in the abelian 
limit is not a severe problem, it should rather be viewed as a technical simplification. 

Corrections to the abelian limit come from two sources: on the one hand, a 
subset of the non-abelian diagrams serves to reconstruct the full one-loop f3 function, 
Po, from its abelian counterpart, (3^; on the other hand, the fact that the event 
shape is sensitive to the details of gluon splitting generates contribution that arc not 
proportional to Pq, since they are not directly related to the running of the coupling. 
The first source of corrections (which must be taken into account also in the simplifed 
massive gluon scheme) is under control: there are strong physical arguments for such 
corrections to be there, and furthermore, at least in principle, they could be explicitly 
included by using existing techniques to isolate gauge-invariant subsets of diagrams 
such as the pinch technique [26]. The second set of corrections is related to the 
splitting of the gluon into two gluons, and these are the corrections included in the 
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soft approximation in the Milan factor. If it were to prove impossible to iron out 
the numerical differences between the two approaches, an analytic computation of 
non-abelian gluon splitting can in principle be performed. 

A more interesting issue from a phenomenological point of view is the study of 
event shape distributions, such as da/dC or daL/dx, with x the energy fraction of the 
detected hadron. Such studies can be performed with the technique outlined here, 
and in fact in [11] the (7^ distribution was studied, showing how the 1/Q correction 
to ai arises from a summation of even power correction to all orders. However, 
having shown that the power correction to the average event shape is determined by 
soft gluon emission, we expect to recover at least qualitatively the results of [12,13], 
namely the characteristic constant shift in the distribution from its perturbative 
estimate by a power suppressed amount. It would be interesting to understand to 
what extent this simple behaviour of the distributions is due to the approximations 
inherent in the dispersive approach. In fact the shift in the distributions is recovered 
in a factorization-based approach [5] only as an approximation of the full answer, 
obtained essentially by neglecting the long-range, wide-angle correlations of the soft 
radiation. A complete description of the power correction to the distribution requires 
in that approach a non-perturbative function rather than a single parameter. The 
phenomenological impact of a more detailed analysis of the power correction along 
these lines is at present an open question. 
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Note After the discrepancy between our results and those of Refs. [12, 13] had 
become evident, an error of a factor of two was found in the computation of the 
Milan factor [27] . Once this error is corrected, there is complete agreement with our 
calculation. Moreover it becomes apparent that ai indeed belongs to the family of 
observables that receive the universal Milan enhancement, there no longer being any 
conflict with Ref. [11]. This should enable experimental investigation of the power 
correction to just as for other shape variables. 
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